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ABSTRACT 

We apply the Batalin-Tyutin Hamiltonian method to the abelian Proca model in 
order to convert a second class constraint system into a first class one systemetically by 
introducing the new field. Then, according to the BFV formalism we obtain that the 
desired resulting Lagrangian preserving standard BRST sysmmetry naturally includes 
the well-known Stiickelberg scalar related to the explicit gauge-breaking effect due to 
the presence of the mass term. Furthermore, we also discuss the nonlocal symmetry 
structure of this model in the context of the nonstandard BRST symmetry. 



*Electronic address: YOONSJ@GSCRL.goldstar.co.kr 
^Electronic address: ywkim@physics. sogang. ac.kr 
■^Electronic address: yjpark@physics. sogang. ac.kr 



1 



1 Introduction 



The Dirac method has been widely used in the Hamiltonian formalism [1] to quan- 
tize second class constraint system, which does not form a closed constraint algebra 
in Poisson brackets. However, since the resulting Dirac brackets are generally field- 
dependent and nonlocal, and have a serious ordering problem between field operators, 
these are under unfavorable circumstances in finding canonically conjugate pairs. On 
the other hand, the quantization of first class constraint system [2,3] has been well 
appreciated in a gauge invariant manner preserving Becci-Rouet-Stora-Tyutin (BRST) 
symmetry [4,5]. If second class constraint system can be converted into first class one 
in an extended phase space, we do not need to define Dirac brackets and then the 
remaining quantization program follows the method of Ref. [2-5]. This procedure has 
been extensively studied by Batalin, Fradkin, and Tyutin [6,7] in the canonical for- 
malism, and applied to various models [8-10] obtaining the Wess-Zumino (WZ) action 
[11,12]. 

Recently, Banerjee [13] has applied the Batalin- Tyutin (BT) Hamiltonian method 
[7] to the second class constraint system of the abelian Chern-Simons (CS) field the- 
ory [14-16], which yields first class constraint algebra in an extended phase space by 
introducing new fields. As a result, he has obtained the new type of an abelian WZ 
action, which cannot be obtained in the usual path-integral framework. Very recently, 
we have quantized several interesting models [17] as well as the nonabelian CS case 
[18], which yield the weakly involutive first class system originating from the second 
class one, by generalizing this BT formalism [7,13]. As shown in these works, the 
nature of the second class constraint algebra originates from the symplectic structure 
of CS term, not due to the local gauge symmetry breaking. Banerjee and Ghosh [19] 
have also considered a massive Maxwell theory, which has the explicit gauge-breaking 
term, in the BT approach. However, all these analyses do not carry out the covariant 
gauge fixing procedure preserving the BRST symmetry. On the other hand, Lavelle 
and McMullan (LM) recently have found that QED exhibits a new nonlocal symmetry 
[20]. Several authors [21,22] have extensively studied following their works. 

In the present paper, we shall apply the BT Hamiltonian method [7] to the Abelian 
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Proca theory revealing the Stiickelberg effect [23]. In section 2, we apply this formalism 
to the abelian Maxwell (Proca) theory in order to systematically convert a second 
class constraint system into a first class one by introducing a new auxilary field p. 
In section 3, we will briefly discuss the special unitary gauge fixing reproducing the 
original second class theory. In section 4, we show that by identifying this unphysical 
new field p with the Stiickelberg scalar we naturally derive the Stiickelberg scalar term 
related to the explicit gauge-breaking mass term through a standard BRST invariant 
gauge fixing procedure according to the Batalin-Fradkin-Vilkovisky (BFV) formalism. 
We also analyse the nonlocal symmetry structure, which exists in QED, of the Proca 
model in the context of the nonstandard BRST symmetry. 



2 The BT Formalism 

Now, we first apply the BT formalism to the abelian massive Maxwell theory in 
four dimensions, whose dynamics are given by 

S = Jd 4 x [~F^F^ + l -m 2 A^A% (1) 

where F^ = d^A u - d u A fM , and = diag(+, -, -, -). 
The canonical momenta of gauge fields are given by 

7T = 0, 

7T, = A l + d t A°. (2) 



Then, Q\ = ttq is a primary constraint [1]. The canonical Hamiltonian is 

H c = J d*x ^ + -F^ + \m\{A*Y + (A*) 2 ) - A Q 2 



(3) 



where fl 2 is the Gauss' law constraint, which comes from the time evolution of Q 1 , 
defined by 

tt 2 = 8% + m 2 A°. (4) 

The time evolution of the Gauss' law constraint generates no more independent con- 
straints. As a result, the full constraints of this model are Qi and Q 2 - Then, they 



consist of the second class constraint algebra as follows 

A ij (x,y) = {n i (x),n j (y)} = -m 2 e ij 6 3 (x-y) = 1,2), (5) 

where we denote e 12 = e 12 = 1. 

We now introduce new auxiliary fields $ l to convert the second class constraint Jl, 
into first class one in the extended phase space, and assume that the Poisson algebra 
of the new fields is given by 

{* i (x),&(y)} = u> i *(x,y), (6) 

where is an antisymmetric matrix. According to the BT method [7], the modified 
constraint in the extended phase space is given by the polynomials of the auxiliary 
fields as follows 

oo 

q^, a», = a + E «i B) ; ~ (*')", (7) 

n=l 

satisfying the boundary condition, Cl^n^, A^, <E> J = 0) = Qj. The first order correction 
term in the infinite series [7] is given by 

tip{x) = jd 3 yX t ,(x,y)^(y), (8) 

and the first class constraint algebra of Cl t requires the condition as follows, 

A tj (x, y)+ J d 3 w d 3 z X ik (x, w)uj k \w, z)X jl (z, y) = 0. (9) 

As was emphasized in Ref. [13,17], there is a natural arbitrariness in choosing and 
Xij from Eq. (6) and Eq. (8), which corresponds to canonical transformation in the 
extended phase space [6,7]. Thus, without any loss of generality, we take the simple 
solutions as 

to ij (x,y) = e ij 8 3 (x-y), 

Xij(x,y) = m5 ij 8 2l (x-y), (10) 

which are compatible with Eq. (9) as it should be. Then, the modified constraint, Cli 
give a strongly first class constraint algebra, 

{Q i (x),Q j (y)} = 0, (11) 
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where 

Qi = a t + m¥ (12) 

are the modified constraints including the auxiliay fields in the extended phase 
space. 

Next, we derive the corresponding involutive Hamiltonian in the extended phase 
space. It is given by the infinite series [7], 

oo 

H = H c + Y^ H (n) ; H {n) ~ (&) n , (13) 

n=l 

satisfying the initial condition, H^^, A^, $ l = 0) = H c . The general solution [7] for 
the involution of H is given by 

= -- f d 3 xd 3 yd 3 z ^(x)u ij (x,y)X^(y,z)G { r 1) (z) (n > 1), (14) 

where the generating functions G^ are given by 

g?> = {n?\H c }, 

Gt ] = {nP,H^}o + {Q^,H^} (n>l), (15) 

where the symbol O in Eq. (15) represents that the Poisson brackets are calculated 
among the original variables, i.e., O = (7r M ,A M ). Here, and X ij are the inverse 
matrices of and X^, respectively. Explicit calculations yield 

G^ = Q2, 

= m 2 diA\ (16) 

which are substituted in (14) to obtain 
if « = J d 3 x 

This is inserted back in Eq. (15) in order to deduce G^p as follows 



m 



(17) 



G? = m$ 2 , 

G ( 2 1] = mdi&Q 1 . (18) 
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Then, we obtain H {2) by substituting Gf ] in Eq. (14) 



= J d 3 x 

Finally, since 



(19) 



= (n > 2), (20) 

we obtain the complete form of the Hamiltonian H after the n = 2 finite truncations 
as follows 

H = H C + H<U +# (2) , (21) 
which, by construction, is strongly involutive, 

{n i ,H} = 0. (22) 

This completes the operatorial conversion of the original second class system with 
Hamiltonian H c and constraints Q{ into the first class with Hamiltonian H and con- 
straints fij. From Eqs. (11) and (22), one can easily see that the original second class 
constraint system is converted into the first class system if one introduces two fields, 
which are conjugated with each other in the extended phase space. Note that the origin 
of second class constraint is due to the explicit gauge symmetry breaking term in the 
action (1). 

Next we consider the partition function of the model in order to present the La- 
grangian corresponding to H in the canonical Hamiltonian formalism. As a result, 
we will unravel the correspondence of the Hamiltonian approach with the well-known 
Stiickelberg's formalism. First, let us identify the new variables as a canonically 
conjugate pair (p, n p ) in the Hamiltonian formalism, i.e., 

$ 1 = (mp, -n p ) (23) 
m 

satisfying Eqs. (6) and (10). Then, the starting phase space partition function is given 
by the Faddeev formula [3,24] as follows 

r 2 
Z = / VA^Vn^VpVTTp [] KQ^diT^det | {fi^Tj | e lS , (24) 
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where 

S = Jd 4 x (tt^ + ix p p - H) , (25) 

with Hamiltonian density H corresponding to Hamiltonian H (21), which is now ex- 
pressed in terms of (p, tt p ) instead of $\ The gauge fixing conditions Tj are chosen so 
that the determinant occurring in the functional measure is nonvanishing. Moreover, 
Tj may be assumed to be independent of the momenta so that these are considered as 
Faddeev-Popov type gauge conditions [24]. 

Before performing the momentum integrations to obtain the partition function in 
the configuration space, it seems appropriate to comment on the involutive Hamil- 
tonian. If we directly use the Hamiltonian (21) following the previous analysis done 
by Banerjee et al. [19], we will finally obtain the non-local action corresponding to 
this Hamiltonian due to the existence of (5 i 7r i ) 2 -term in the action when we carry out 
the functional integration over n p later. Furthermore, if we use this Hamiltonian, we 
can not also naturally generate the first class Gauss' law constraint Vt 2 from the time 
evolution of the primary constraint Qi, which is the first class. Therefore, in order to 
avoid these serious problems, we use the equivalent first class Hamiltonian without any 
loss of generality, which only differs from the involutive Hamiltonian (21) by adding a 
term proportional to the first class constraint Q 2 as follows 

H' = H+^n 2 . (26) 

Then, we have the natural first constraint system such that 

{n 1 ,H'} = n 2 , {n 2 ,H'} = o. (27) 

Note that when we act this modified Hamiltonian (26) on physical states, the difference 
is trivial because such states are annihilated by the first class constraints. Similarly, 
the equations of motion for observable (i.e. gauge invariant variables) will also be 
unaffected by this difference since Vl 2 can be regarded as the generator of the gauge 
transformations. 
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3 The Original Unitary Gauge Fixing 



Now, we consider the following effective phase space partition function 
r 2 

i,3=l 

S' = Jd A x (n^ + TTpp-H'). (28) 



The trivial 7r integral is performed by exploiting the delta function 5(£li) = <5(7r +m 2 p) 
in (28). On the other hand, the other delta function 5(f2 2 ) = 8(d l iTi + m 2 A° + tt p ) can 
be expressed by its Fourier transform with Fourier variable £ as follows 

5(n 2 ) = Jv^e-if** (29) 

Making a change of variable A — > A + £, we obtain the action 

S = Jd 4 x - m 2 p(> + + 7r p p - |tt? - ^FijF ij + \m 2 (Ay + ^m 2 A^ 
+ - m 2 ^ - ^tt 2 + im^pdV " f - \m 2 e\ , (30) 

where the corresponding measure is given by 

[Dp] = V^V^iDA^ViipV pY\{5[T j(A° + £, A\ tt,, p, n p )]}det \ {Q h T 3 } \ . (31) 



Performing the Gaussian integral over 7Tj, this yields the action as follows 
S u = J d A x 



\F^ U + lrn 2 A»A» + n p (p-S-^- 2 



- m 2 p(A° + £)- n^d^p + \m 2 d i pd i p - -m 2 £ 2 



(32) 



Now, we choose the unitary gauge as follows 

r i = (P, n p ). (33) 

Note that this gauge fixing is consistent because when we take the gauge fixing condition 
p — 0, another condition % p = is naturally generated from the time evolution of p, i.e., 
p — {p, H u } = —j^np = 0, where the Hamiltonian H u corresponds to the intermediate 
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action S u . However, we can not smoothly choose the massless limit in this unitary 
gauge because p tends to infinity for this limit. In this gauge, we get 

{n t (x),T J (y)} = e tJ 5 3 (x-y). (34) 

Then, we easily recover the original system. Therefore, we can interpret the original 
system (1) as a gauge-fixed version of the extended gauge system (13), (21), and (26). 



4 The BFV-BRST Gauge Fixing 

In this section, we first briefly recapitulate the BFV formalism [2,3] which is appli- 
cable for the general theories with first-class constraints. For simplicity, this formalism 
is restricted to a finite number of phase space variables. This makes the discussion 
simpler and conclusions more apparent. 

First of all, consider a phase space of canonical variables q l , pi (i = 1, 2, • • •, n) in 
terms of which the canonical Hamiltonian H c (q l ,pi) and the constraints Vt a (q l ,pi) « 
(a — 1, 2, • • •, m) are given. We assume that the constraints satisfy the constraint 
algebra [2,3] 

[n a ,n b ] = tn c u c ab , 

[H c ,Q a ] = iQ b V a b , (35) 

where the structure coefficients U c ah and V% are functions of the canonical variables. We 
also assume that the constraints are irreducible, which means that locally there exists 
an invertible change of variables such that Q a can be identified with the m-unphysical 
momenta. 

In order to single out the physical variables, we can introduce the additional condi- 
tions Q a (q\pi) ~ with det[$ a , Q b ] ^ at least in the vicinity of the constraint surface 
$ a and Vl a ~ 0. Then, $ a play the roles of gauge- fixing functions. That is to say, 
from the condition of the time stability of the constraints, there exists a family of 
phase space trajectories. By selecting one of these trajectories through the conditions 
of $ a ~ o, we can get the 2{n — m) dimensional physical phase space noted by q* and 
p* [1-3]. And then, Q a (q\pi) can be identified with the m-unphysical coordinates. 
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The described dynamical system 

Z = J[dq l d Pl ] 5(tt a )5(<5> b )det | [<5> b ,Q a ] | e ^ dx{p ^ Hc) (36) 

is completely equivalent to an effective quantum theory only depending on the physical 
canonical variables q,p of the physical phase space [2,3]. Note that the constraints 
Q a ~ and $ a ~ together with the Hamilton equations may be obtained from a 
action 

S = jdt ( Pl q l -H c - N a il a - B a <$> a ) , (37) 

where N a and B a are Lagrange multiplier fields canonically conjugated to each other, 
obeying the commutation relations 

[N a ,B b ] = i5 a b , (38) 

and the gauge-fixing conditions contain A a in the following general form 

$ a = N a + x a (q\PhN a ) , (39) 

where \ a are arbitrary functions. Furthermore, we can see that the Lagrange multi- 
plier N a become dynamically active, and B a serve as their conjugate momenta. This 
consideration naturally leads to the canonical formalism in an extended phase space. 

In order to make the equivalence to the initial theory with constraints in the reduced 
phase space, we may introduce two sets of canonically conjugate, anticommuting ghost 
coordinates and momenta C a , V a and V a , C a such that 

[C a ,V b ] = [V a ,C b ] = i5 a b . (40) 

The quantum theory is defined by the extended phase space functional integral 

Z* = J [dfi]e iS * , (41) 

where the action is now 

S* = Jdt{ Pi q l + B a N a + V a C a + C a V a - H m + i[Q,*]}, (42) 
and [d/j] is the Liouville measure, i.e., 

[dfi] = [dq i d Pi dN a dB a dC a dV a dV a dC a \, (43) 
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on the constraint phase space. Here, the BRST-charge Q and the fermionic gauge-fixing 
function \I> are defined by 

Q = c a n a - l -c b c c u%v a + v a B a , 

V = C aX a + V a N a , (44) 

respectively. H m is the BRST invariant Hamiltonian, called the minimal Hamiltonian, 

H m = H c + V a V b a C b . (45) 

Now, in order to derive a BRST gauge-fixed covariant action for the abelian Proca 
model considered in the previous section, let us introduce the ghosts and anti-ghosts 
along with auxiliary fields (C\ Vi), (V\ C,), (N\Bi), where i — 1, 2, according to 
the above BFV formalism in the extended phase space. The nilpotent BRST-charge 
Q, the fermionic gauge-fixing function and the minimal Hamiltonian H rn have the 
following concrete forms 

Q = Jd 3 x[ C% + V i B i ], 
= \<$x\C lX l + ViN*], 
H m = H' — Jd 3 x \p 2 C\ (46) 

where 

X l = A\ X ^ = d^ + ^B 2 , (47) 

as gauge fixing functions, and a is an arbitrary parameter. Note that the form of H m 
is simpler than that in Ref. [8] due to our improved Hamiltonian (26) proposed in our 
previous works [17]. 

The BRST-charge Q, the fermionic gauge-fixing function and the minimal Hamil- 
tonian H m satisfy the following relations, 

i[Q,H m ] = 0, 
Q 2 = [Q,Q] = 0, 

[[%Q],Q] = 0, (48) 
where they are the conditions of physical subspace being imposed. 
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Then, the effective action is 

S eff = Jd^i n A° + 7r t A l + n p p + B 2 N 2 + V t C + C 2 V 2 ] - H total , (49) 

where H tota i = H m — i[C t ), x I / ]. Note that we could suppress the term / d 4 x(BiN l + 
C{P V ) = —i[Q,f d A xC\N l ] in the Legendre transformation by replacing x 1 with x 1 + 
N 1 . 



4.1 The Standard Local Gauge Fixing 

The fields B 1 , N 1 , Ci, V 1 , Pi, C 1 , and A are eliminated, and integration of 7r gives 
the delta functional by using of Gaussian integration. Then we obtain the generating 
functional as follows 

Z = J[Dp]exp[tS eff ], (50) 

where the effective action is 

1, 1 
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S eff = J d A x [ mA l + ix p p + BN + VC + CV - -{n 1 ) 2 - 

--F^ - l -m 2 (A*) 2 - rrfpdiA* + l -m 2 d t pd l p 

+N(d l n i + u p ) + B{d l A l + ^aB) + diCd'C + VV], (51) 

and the Liouville measure of the extended phase space is given by 

[Dp] = [Dm VA i Vp Vn p VB VN VV VC VC VV]. (52) 

Here we have redefined N 2 = N, B 2 = B, C 2 =C,C 2 = C, V 2 = V, and V 2 = V. 
Performing the integrations of ir\ ti p , V, and V, and identifying with iV = —A , we 
get the following covariant effective action 

S eff = J d A x \--F pu F^ + \m\A p + d,p) 2 - A»d,B + l -a{B) 2 - d p Cd»C ], (53) 

which is invariant under the BRST transformation 

5 B A p = —Xd p C, S B p = \C, 
S B C = 0, 5 B C = -XB, 5 B B = 0, (54) 
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where A is a constant Grassmann parameter, and the corresponding final measure is 
given by 

[Dp] = [VA» Vp VB VC VC\. (55) 

Therefore, in Eq. (53) we see that the auxiliary BF field p is exactly the well-known 
Stiickelberg scalar [23]. Note that in this gauge, we can smoothly choose the massless 
limit because the gauge fixing conditions (47) are independent of mass m 2 , and then 
obtain the well-known QED result. On the other hand, this BRST symmetry gives a 
conserved current as 

Jbh = F^cTC + m\A^ + 8^p)C + Bd^C, (56) 
through Noether's theorem. 



4.2 The Nonstandard Nonlocal Gauge Fixing 

Consider the BFV formalism in the previous section up to the point, where the in- 
tegration over the momentum n p was performed and the following action was obtained. 

S eff = Jd 4 x + Kn\A, + d,p) 2 - A»d,B + l -aB 2 

-d£d l C + VC + CV + VV J . (57) 

This action is invariant under the BRST transformation, which have the form 

5 B A = -XV, S B A t = -Xd t C, 5 B p=XC, 
5 B C = 0, 5 B C = -XB, S B B = 0, 

5 B V = 0, 5 B V = -X[-d l F 0l + m 2 {p + A )]. (58) 

Now, if we perform the integration over the ghost fields instead of their momenta, 
we can also find the nonlocal symmetry structure in the Proca theory as well as QED 
[20]. First, performing the integration over the ghost field C, we get the following delta 
function 

Sidid'C-V') = det^did^C-^—V). (59) 
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Next, performing the integration over C, we get non-local ghost Lagrangian 

S gh = Jd'x^^V-VV]. (60) 

Notice that the appearance of the nonlocal term in the ghost action has a result of this 
unusual integration. However, this form can be also simply obtained by the change of 
variables 

c ^h*'~ c ^ v ' (61) 

in the Eq. (53). Under these replacements, we have the nonlocal BRST charge, that 
we call Q', being 




Then, the effective action is BRST invariant under the following transformations as 

5' B V = 0, 5' B V = -XB, 5' B B = 0. (63) 

This nonlocal BRST symmetry yields a conserved current through Noether's theorem 
as follows 

J' B , = + ™ 2 <A. + M^f* + B%±P. (64) 

Note that performing the change of variable (61), these nonlocal symmetry and con- 
served current turn into just the original local theory (56). 

In conclusion, we have applied our improved Batalin-Tyutin method [17], which 
systematically converts the second class system into the first class one by introducing 
the new auxiliary fields, to the Proca theory. According to the BFV formalism with 
the efficient first class Hamiltonian through BT analysis, we have shown that the 
resulting Proca Lagrangian preserving standard BRST symmetry naturally includes 
the well-known Stiikelberg scalar needed for the anomaly free theory by identifying 
this scalar with one of auxiliary fields. Furthermore, we have shown that the nonlocal 
symmetry structure recently proposed in QED also exists in the Proca model through 
the nonstandard BRST gauge-fixing procedure. 
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